Direct Measurement of the Flow Field Around Swimming Microorganisms 
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Swimming microorganisms create flows that influence their mutual interactions and modify the 
rheology of their suspensions. While extensively studied theoretically, these flows have not been 
measured in detail around any freely- swimming microorganism. We report such measurements for 
the microphytes Volvox carteri and Chlamydomonas reinhardtii. The minute (~ 0.3%) density 
excess of V. carteri over water leads to a strongly dominant Stokeslet contribution, with the widely- 
assumed stresslet flow only a correction to the subleading source dipole term. This implies that 
suspensions of V. carteri have features similar to suspensions of sediment ing particles. The flow in 
the region around C reinhardtii where significant hydrodynamic interaction is likely to occur differs 
qualitatively from a "puller" stresslet, and can be described by a simple three- Stokeslet model. 

PACS numbers: 87.17. Jj,87.16.Qp,47.63.Gd 



Aided by advances in imaging techniques that allow de- 
tailed studies of the rotating flagella of bacteria [l| and 
the undulating flagella of spermatozoa [2] and algae 
there is now a general consensus on how mechanical mo- 
tions of microorganism appendages generate propulsive 
forces in a viscous fluid [4] . No such consensus exists yet 
on the origins of collective behavior ^ , transport @, Q 
and rheological properties of suspensio ns I8| , and the in- 
teraction of organisms with surfaces [9|, LLOj . As hydrody- 
namics surely plays a key role in these effects, a detailed 
knowledge of the flow field around freely swimming mi- 
croorganisms is needed, both in the near-field and far 
away. Here we present the first such measurements. 

The linearity of the Stokes equations implies that the 
far-field flow around a microorganism can be expressed 
as a superposition of singularity solutions [11], with the 
slowest decaying mode dominating sufficiently far away. 
Theories of fluid-mediated interactions and collective 
behavior typically assume neutrally buoyant swimmers 
which exert no net force on the fluid. The thrust T of 
their flagella and the viscous drag on their body are dis- 
placed a distance d apart (often comparable to the cell 
radius R) , and balance to give the far-field flow of a force 
dipole, or stresslet [l2j |. which decays with distance r as 
Td/r]r 2 , where n is the fluid's viscosity. The contribution 
from a suspension of such stresslets to the fluid stress ten- 
sor is central to some of the most promising approaches 
to collective behavior of microorganisms [13|- 

The force-free idealization of swimmers requires pre- 
cise density-matching [9] not generally realized in na- 
ture. To appreciate the striking effects of gravity, one 
need only consider the buoyancy-driven plumes of bio- 
convection [l4|. Models of this instability express the 
contribution of cells to the Navier- Stokes equations as 
a sum of force monopoles (Stokeslets), coarse-grained as 
a body force proportional to the cell concentration and 
gravitational force F g per cell [14]. As the flow around a 
Stokeslet decays as Fg/nr, it is clear, if not appreciated 



previously, that there is a distance A ~ Td/F g at which 
the nearby stresslet contribution crosses over to the dis- 
tant Stokeslet regime. This is one of several crossover 
lengths relevant to swimmers; for ciliates, unsteady ef- 
fects become important on scales smaller than the vis- 
cous penetration depth [15[. For a given organism, the 
relevance of the length A to a particular physical situ- 
ation depends on the cell concentration and the observ- 
able of interest. At low concentrations the Stokeslet form 
suffices, but the near field is relevant to cell-cell interac- 
tions, especially in concentrated suspensions [16| and to 
tracer dynamics [7[ . The notion of near field requires dis- 
tinguishing between distances r satisfying i? « r « A, 
where a stresslet description may hold, and r ~ R where 
the multipole contributions may not be well-ordered and 
the flow topology can differ from that of a stresslet. 

A synthesis of tracking microscopy and fluid velocime- 
try is used here to quantify the competing force singu- 
larities and the near-field flow topology for the unicel- 
lular biflagellate green alga Chlamydomonas reinhardtii 
171 (R ~ 5 /im) and its larger relative Volvox carteri 
18|, a spherical alga (R ~ 200 /im) which swims by the 
action of ~ 10 3 Chlamydomonas -like cells on its surface. 
For Volvox our most significant finding is that the flow 
field is strongly dominated by its Stokeslet component, 
despite a density excess of a mere ~ 0.3%, much smaller 
than that of common unicellular organisms (~ 5 — 10%). 
Moreover, the high symmetry of Volvox results in a lead- 
ing near-field correction in the form of a source doublet, 
and a smaller stresslet. The flow around Chlamydomonas 
is compatible with a simple "puller" stresslet only at dis- 
tances > 7R, where the fluid velocity is < 1% of the 
swimming speed; closer to the cell, the flow topology re- 
flects the finite separation of the flagellar and body forces. 

V. carteri f. naqariensis (strain EVE) was grown ax- 
enically in SVM [19] with sterile air bubbling, whereas 
C. reinhardtii (strain UTEX 89) was grown axenically in 
TAP medium [17] on an orbital shaker, both in a diurnal 



2 




FIG. 1: (color online). Flow field of a freely swimming V. 
carteri in the laboratory frame. (a,b) Magnitude and stream- 
lines of u and its fitted approximation Uf it respectively, (c) 
Relative error of the fit: e = ||u — Uf it ||/||u||. (a) and (b) have 
the same colorbar, different from (c). g indicates gravity. 



growth chamber with 16 h in artificial cool daylight (~ 
4000 lux) at 28° C, and 8 h in the dark at 26° C. The large 
difference in organism size between Volvox and Chlamy- 
domonas required two distinct methods to measure the 
flows they create [20]. A CCD camera (Pike, Allied Vi- 
sion Technologies) mounted on a continuously-focusable 
microscope (Infinivar, Infinity Optics) and connected to 
a vertical motorized XY stage (Thorlabs) followed indi- 
vidual Volvox colonies as they swam upwards [2l| in a 
straight line along the central axis of a 5 x 5 x 50 mm 
sample chamber filled with SVM at 21 ± 1° C. The stage 
was controlled by a custom Lab View routine. The fluid 
was seeded at volume fraction 10 ~ 5 with 1 jam nile-red 
polystyrene microspheres (Invitrogen) illuminated by a 
vertical ~ 500 jam thick laser sheet (A = 532 nm). Volvox 
is phototactic [22] at this wavelength, and at the inten- 
sities used here it swims smoothly along the laser sheet. 
We recorded the flow field of 19 different colonies at 30 fps 
for ~ 2 — 3 min each. The measured flow field v was ob- 
tained by particle image velocimetry (Dantec Dynamics). 
Background flows in the chamber were < 10 /im/s. 

We observed a dilute suspension (~ 3 x 10 6 cells/cm 3 ) 
of Chlamydomonas in TAP on a Nikon inverted micro- 



scope at 40 x (NA 0.6) by exciting their chlorophyll aut- 
ofluorescence with a laser (635 nm, ~ 60 mW), which 
also excited 1.6 pm fluorescent polystyrene microspheres 
(FS04F, Bangs Labs) used as tracers. Cylindrical poly- 
dimethylsiloxane sample chambers (5 mm radius, 0.4 mm 
height) were prepared, pacified, and filled following 
Experiments were performed at 21 ± 1° C, with the laser 
providing the only light source. We focused on a plane 
150 /im inside the chamber to minimize surface effects, 
and recorded movies at 250 fps (Fastcam SA3, Photron). 
Movies were analysed with standard algorithms to track 
cells and tracers. For each cell swimming along the focal 
plane for more than 1 s (~ 10 body lengths), we collected 
the instantaneous velocity of all tracers at r < 14 R, nor- 
malized by the swimmer's speed. The resulting 3.3 x 10 6 
velocity vectors were binned into a 2.5 pm square grid 
(shown in Fig. [4] below), and the mean of the well- 
resolved Gaussian in each bin was used for the flow field. 

In both experiments Uo indicates the swimmer veloc- 
ity, while u(r) and v(r) = u(r) — Uo are the velocity field 
in the laboratory and comoving frames respectively. 

A typical experimental flow field around Volvox is 
shown in Fig. Ufa). We fit these fields to a superpo- 
sition of a uniform background velocity (Uo), a Stokeslet 
(St), a stresslet (str) and a source doublet (sd): 
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where I is the unit tensor, y is the upward vertical unit 
vector, r = r/r, and r is measured from the center of 
the organism (x c , y c ). The orientation of all multipoles is 
fixed to be along the vertical, and we are left with six pa- 
rameters: (/To, Ast, A str , A s d , x c , y c ). The fits, obtained 
by minimizing the integrated squared difference between 
the model and the experimental flows, describe remark- 
ably well the experimental flow, almost down to the sur- 
face of the organisms [see Fig. [Hb,c)]. Typical values 
for the parameters are Uo ~ 10 2 /im/s, Ast ~ 10 4 /im 2 /s, 
A str ~ 10 6 /im 3 /s (indicating a pusher- type stresslet), 
A s d ~ 10 9 /im 4 /s, with the actual magnitude depend- 
ing on the colony radius R. From the Stokeslet com- 
ponent, we can calculate the average colony density as 
Ap = 6rjAst/gR 3 , where rj = 10 _3 Pas and g is the grav- 
itational acceleration. The dependence of both Ap and 
Uo on R (Fig. [2j) compares well with previously published 
data [21] obtained by different means, thereby validat- 
ing the measurements and analysis procedures. Remov- 
ing the Stokeslet contribution from the experimental flow 
field [Fig. [3ja)] reveals that the near field is dominated by 
the source doublet component, with the stresslet respon- 
sible only for a slight forward-backward asymmetry [Fig. 
|3fb,c)]. The orientation of the source doublet is opposite 
to that around a translating solid sphere, and is compat- 
ible with a model that assigns a constant force density to 



3 





100 200 300 400 500 600 700 
ft (Mm) 



FIG. 2: (color online). The dependence of the excess den- 
sity Ap and swimming speed Uo (inset) of V. carteri colonies 
(red stars) on their radius R are compatible with previous 
measurements (blue circles) [21], and also [251 ]. 



the colony surface |23], as well as a particular case of the 
"squirmer" model [24]. Average values of the parameters 
in Eq. [T] show that the crossover distance between the 
source doublet and stresslet (~ 3 R) is beyond that at 
which the Stokeslet becomes the leading component of 
the flow (~ 1.5 R). This peculiar ordering of multipoles 
results from the high degree of anterior-posterior sym- 
metry of Volvox [23j, and highlights the influence that a 
swimmer's body plan can have on its flow field. 

The two ~ 12 fim long flagella of Chlamydomonas 
beat mostly in a synchronous breast stroke at ~ 50 Hz 
0|, pulling the cell body through the fluid at speeds 
Uo ~ 100/im/s. Despite its ~ 5% density excess over 
water [25[ , the gravitational Stokeslet of Chlamydomonas 
only becomes dominant at distances A > 35 R, as its ra- 
tio of Uo to the sedimentation speed (F g /67rr]R) is much 
larger than for Volvox. Therefore this swimmer has often 
been modelled as a puller stresslet [l4|. A slight three- 
dimensionality of the beating causes the cell to spin about 
its swimming direction at ~ 2 Hz, so ensemble- averaged 
measurements of the kind presented here average out az- 
imuthal asymmetries in the flow field. Figure [4f a) shows 
that for r > 7R the measured flow topology begins to 
resemble a puller stresslet, yet flow speeds at such dis- 
tances are already < 1 /im/s. Closer to the organism, 
the field becomes more complex. It includes side vortices 
and a flow in front of the cell body that is along the direc- 
tion of motion, towards a stagnation point. The velocity 
field can be accurately captured, even in the near-field, 
by modeling the flow created by the pulled cell body as 
a Stokeslet, distributing the thrust among two Stokeslets 
located at the approximate positions of the two flagella, 
and averaging the flow over one rotation about Uo- The 
streamlines of this simple extension to the force dipole 
model [see Fig. IH[b)], as well as the decay of ||u(r)|| with 




FIG. 3: (color online). Near field around V. carteri. (a) 
Magnitude, vector fields and streamlines of u after subtract- 
ing the fitted Stokeslet. Colorbar as in Fig. [TJa). (b) ||u|| 
along a horizontal section through the center of the organism. 
The average Stokeslet (green dashed line) follows the decay 
of the experimental flow (blue circles) averaged over 19 differ- 
ent colonies (black dots). Deviations from a pure monopole 
appear from < 5R, and can be captured adding a source dou- 
blet and a stresslet (red solid line), (c) Vertical section of the 
flow field u from the experiment in (a) through the center 
of the colony. The stresslet component is responsible for the 
forward-backward asymmetry. Symbols as in (b). 



distance [see Fig. H(c)], are very similar to those mea- 
sured. Including no-slip boundary conditions on the cell 
body [26[ has little effect on the velocity field as the cell- 
drag Stokeslet nearly produces the appropriate velocity 
field on the cell surface. 

These flow field measurements around freely-swimming 
microorganisms provide the basis for a deeper under- 
standing of a number of issues in biological fluid dynam- 
ics, including the interactions of microorganisms with 
surfaces, with each other, and the rheology of suspen- 
sions. For example, it was recently discovered that Volvox 
colonies can form hydrodynamic bound states whose 
properties are quantitatively described by a model of in- 
teracting Stokeslets near a no- slip wall Hl|. 

The near 

complete dominance of the flow field around Volvox by 
the Stokeslet term found here provides ex post facto justi- 
fication for the neglect of higher moments. Perhaps more 
importantly this result shows that in terms of interparti- 
cle hydrodynamic couplings a suspension of Volvox is like 
a sedimenting suspension [27J , except that the velocity 
of each colony is the sum of a self-propelled contribution 
and mutual advection in the flow field of other spheres. 
Elsewhere we illustrate this correspondence in detail [HI . 

The correspondence between the measured time- and 
azimuthally-averaged flow field of Chlamydomonas and 
the three- Stokeslet model illustrates how well such a sim- 
plification captures the complex flow topology, lending 
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FIG. 4: (color online). Time- and azimuthally-averaged flow field of C reinhardtii. (a) Streamlines (red) computed from 
velocity vectors (blue). The spiraling near elliptic points is an artifact of the direct integration a noisy experimental velocity 
field. A color scheme indicates flow speed magnitudes, (b) Streamlines of the azimuthally averaged flow of the three- Stokeslet 
model: flagellar thrust is distributed among two Stokeslets placed (not fitted) at the approximate flagellar position (lateral 
green arrows), whose sum balances drag on the cell body (central red arrow), (c) Decay of ||u(r)|| for the three directions 
indicated by separate colors in the inset, compared to results from the three-Stokeslet model (dashed). 



support to this approximation in modelling ciliary inter- 
actions [29(| . Our results indicate that the simple puller- 
type description for Chlamydomonas is only valid at dis- 
tances > 7 R, where the flow field is already < 1% of 
Uq. We then expect interactions with other swimmers, 
boundaries or tracers, to be influenced mostly by the flow 
structure at shorter separations, where the full time de- 
pendence of the flow may be important [30] . We are cur- 
rently investigating whether similar conclusions hold for 
the flow field around bacteria, the prototypical "pusher" 
microorganisms. 
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support from the EPSRC, the BBSRC, the Marie-Curie 
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